Is state tomography an unambiguous test of quantum entanglement? 
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We provide an alternative interpretation of recently published experimental results that were rep- 
resented as demonstrating entanglement between two macroscopic quantum Josephson oscillators. 
We model the experimental system using the well-established classical equivalent circuit of a resis- 
tively and capacitively shunted junction. Simulation results are used to generate the corresponding 
density matrix, which is strikingly similar to the previously published matrix that has been declared 
to be an unambiguous demonstration of quantum entanglement. Since our data are generated by a 
classical model, we therefore submit that state tomography cannot be used to determine absolutely 
whether or not quantum entanglement has taken place. Analytical arguments are given for why the 
classical analysis provides an adequate explanation of the experimental results. 
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The possibility that the supposed quantum behavior of 
a physical system might be interpreted within a classical 
framevi^ork has been of interest for decades. For example, 
regarding Peicrls transitions in quasi one-dimensional 
metals, it was suggested that depinning occurring due to 
quantum tunnelling could explain experimental results 
[1], but it was also found that the same results could 
likewise be explained by modelling the metal classically 
with internal degrees of freedom [2]. Similar observa- 
tions of dual interpretations were made in regards to the 
Pancharatnam-Berry's phase in optical systems, where 
an extra geometric phase [3] arises from adiabatic cyclic 
evolution of a physical system. In this case the quan- 
tum interpretation of the phase [4] was re-interpreted 
by several authors [5,6]. With respect to superconduct- 
ing circuits that include Josephson junctions there are 
two conceptual approaches to interpreting experimen- 
tal outcomes: The classical Resistively and Capacitively 
Shunted Junction (RCSJ) model [7,8], and the quantized 
model put forward by A. J. Leggett and co-workers [9]. 
According to the latter, under appropriate conditions of 
temperature and biasing, the phase difference across a 
Josephson junction would behave as a macroscopic quan- 
tum variable [10,11,12]. In contrast to the classical RCSJ 
model, the quantized model thereby assumes intrinsic 
discrete energy states in the Josephson potential well, 
and escape from the well is possible by quantum tunnel- 
ing of the phase variable through the potential barrier. 

The first observation of Josephson quantum behavior 
for a system operated in washboard potential wells [7,8] 
was reported in 1985 [11]. Thus began two decades of 
investigations of Josephson systems addressing the in- 
triguing possibility of macroscopic quantum behavior as 
it would appear with respect to microwave-induced res- 
onant tunneling, Rabi-oscillations, Ramsey- fringes, and 
spin-echo [12]. The essential concept in explaining these 
experimental observations was the spacing between the 
discrete intrinsic energy levels in the shallow well in com- 
bination with the energy carried into the system by ap- 
plied microwaves. The possibility of a dual interpretation 
of these observations has been pursued [13] by a system- 
atic reconsideration of the experiments from the perspec- 
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FIG. 1: Equivalent circuit of capacitively coupled qubits. 
Each superconducting loop of inductance L contains a Joseph- 
son junction modeled by the RCSJ subcircuit. (b) Potential 
energy (in units of /c$o/27r) of a superconducting loop inter- 
rupted by a Josephson junction as a function of the junction 
phase; the assumed normalized flux bias is 0.6941. The "vir- 
tual particle" representing the system is shown sitting at the 
bottom of its shallow well. 



five of the RCSJ model. We have found that the clas- 
sical description gives a good agreement with reported 
experiments on resonant escape, Rabi oscillations, Ram- 
sey fringes, and spin echo in Josephson systems, all to an 
impressive degree of accuracy, given the simplicity of the 
RCSJ circuit model. Recently [14] quantum mechanical 
entanglement was reported for a system of two weakly 
coupled superconducting qubits. The significance of this 
report is that entanglement is exclusively a quantum phe- 
nomenon. It was claimed that "A full and unambiguous 
test of entanglement comes from state tomography." It 
is that assertion we address in this paper. More partic- 
ularly, we ask the question: From the perspective of the 
classical model, what would tomography say? 

The experimental configuration of two coupled qubits 
is modeled with the equivalent circuit depicted in 
Fig. 1(a). Each junction is characterized by a critical 
current /c, resistance R, and capacitance C; each loop 
has an inductance L. The two loops, characterized by 
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the Josephson phase variables ipi and (p2, arc coupled 
through a capacitance C'x ■ The externally applied fluxes 
through loops 1 and 2 arc represented by and ^2x- 
As discussed in [15], in dimensionless form the governing 
equations of this system turn out to be: 

If a + aipa + sin (fia COS ipb = -/3£^ [fa + '^nMa] (1) 

g~^'fib + a'fib + sin ifij, COS (pa = -Pl^ [fb + '^T^Mb] , (2) 

with overdots denoting derivatives in dimensionless time 
r = ujjt. ifa = {(fi + <P2)/2 and ifb = {^fi - (p2)/2 
are the transformed phase variables that are convenient 
for the analysis. The junction plasma frequency is 
ujj = y/2eIc/hC; a = l/ujCR and Pl = 27rL/c/$o, 
$0 = h/2e being the flux quantum. = 1 + 2"^^ 

with 7^ = Cx/C representing the mutual coupling co- 
efflcient. Adi^ = ^ix/^o is the total normalized ap- 
plied flux through loop i; Ma = ^{Mi^ + M2x) and 
Mb = ^{Mix — M2x) are corresponding transformed nor- 
malized magnetic fields. As in [15], parameter values 
were set at a = 5 x 10-^ Pl = 2.841, g ^ 0.9954, 
Ic = I.IM, C = 1.3pF, and wj^ = 0.02ns. Both 
loops are biased with a dc flux Mi^ = 0.6941 (resonance 
« 5.1GHz) and with superimposed microwave (MW) 
pulses as shown in Fig. 3(A) of Rcf. [14]. 

Our previous work [15] outlined the dynamical modes 
exhibited by this system. One of these modes is con- 
sistent with the experimental observations and can be 
obtained in the linear limit of oscillation amplitude. Fol- 
lowing Ref. [15], we therefore take Mb = 0, a w 0, and 
for small amplitude oscillations of ipi and ip2, we write 
fa = fo+ipa, where ipo is a constant and \ipa\ <C 1. 
Similarly we assume \(pb\ <C 1. Inserting this ansatz with 
a = (for simplicity) gives the linear equations 

/3Lsin(^o + </'o + 27rMa = (3) 

4>a + {cOS(po+P2^)i'a = (4) 

+ (cos(/?o = 0- (5) 

The first of these equations provides the average phase 
ifio for each of the two loops given the external DC 
magnetic field. The next two equations have solutions 
i/ja = Asiii{ujat + Oa) and ipb = Bs\n.{ujbt + Ob), where 
the four parameters, A,B, 9a and 9b are given by ini- 
tial conditions. The two resonance frequencies are = 
cos</5o + Pl^ ^^"^ '^b = \/g^^a- The Josephson phases ipi 
of the two loops therefore evolve according to 

^Pi-VO = V'a ± V>b (6) 

= + B^± 2AB cos2{u)dt + 9d) sin(a;^i -|- 9i) , 

where i = 1,2, and and "— " apply to i = 1 and 
i = 2, respectively. The frequencies and phases are 
= |(wa + Wfe), LOd = |(wa - Wft), and 9d = ^{9a- 9b). 
This expression directly gives the experimentally [14] and 
theoretically [15] observed slow envelope modulation of 
the two coupled oscillators, and calculating the modu- 
lation frequency from the experimentally provided sys- 
tem parameters yields ftj = 2cod ~ 0.001455, which is in 



very good agreement with observations [14,15]. The two 
phases Oi and 92 are also modulated with ujd such that 

sin^i = Acos{u!dt + 6a,) ± B cos{ojdt — 9b) (7) 
cos^j = —Asm{u)dt + 9a)±Bsm{uidt — 9b)- (8) 

We now discuss how these simple linear oscillations 
are related to the Bloch-vectors that have been used in 
the literature to illustrate the system behavior. In anal- 
ogy with the quantum mechanical picture, we define a 
Bloch-sphere centered in a Cartesian coordinate system 
with a horizontal a;y-plane and a vertical z-axis. The cor- 
responding Bloch vector of an oscillator is given by the 
phase and amplitude of oscillation such that the phase 
of oscillation provides the direction in the xy-plane and 
the amplitude (or, equivalently, the energy) provides the 
2;-coordinate. The definition of the vertical axis is given 
by the switching probabilities when the probe pulse is 
applied: z = —1 (vertical down) is a state with near- 
zero switching probability; z = 1 (vertical up) is a state 
with near-unity switching probability; and z = (hor- 
izontal) are states with intermediate switching proba- 
bility. In the present system of two coupled Josephson 
loops, the z = states (ideally 50% switching) are gener- 
ated by initially applying a tt-MW pulse of half a Rabi- 
period duration (notice that Rabi-type oscillations are 
already characterized in the classical system [13]) to one 
of the loops, then allowing the two coupled oscillators 
to exchange energy during a time of free evolution until 
their energy contents, and therefore switching probabil- 
ities, coincide (see Refs. [14,15]). Notice that the ex- 
periments can only measure the vertical axis z by the 
switching probability. Thus, in order to experimentally 
gain any phase information from an oscillator (i.e., the 
■Ty-coordinatc of the Bloch vector in the xy-plane), one 
must "rotate" the Bloch vector ^ around horizontal axes 
of the Bloch-sphere before conducting the measurement 
[14]. This is done by applying additional J-MW pulses of 
a duration that is about half of the initial tt-MW pulse. 
The phase-difference between such rotation pulse and the 
perturbed oscillator determines if the oscillator energy is 
amplified {z increases) or attenuated {z decreases). Ap- 
plying pulses of different phases can therefore illuminate 
the phase (the xy coordinates) of a Bloch-vector by sub- 
sequent switching (z) measurements. A pair of two such 
orthogonal j pulses are denoted X and Y. This is the 
classical analog to state tomography. 

We conduct state tomography by simultaneous switch- 
ing measurements of the two loops after having applied 
rotation pulses (/, X, Y). We acquire the probabilities in 
the form, where, e.g., PxiYq is the probability for simul- 
taneously measuring a switch in loop 1 and a no-switch in 
loop 2 after having applied an X-rotation to loop 1 and 
a y-rotation to loop 2. The notation for not applying an 
X 01 Y rotation is /. Following the quantum mechanical 
treatment, these measured probabilities can be fed into 
expressions for the Hermitian density matrix {pjk}4x4 
in, e.g., the following way: pn = Pioio, P22 = Pioh, 



3 



P33 = Phia, P44 = Phh, and 
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Our only deviation from the experimental work [14] is 
that we use the direct expressions above instead of a 
least-square fit. 

Simulations were conducted as follows. Each loop was 
initially at rest in its shallow potential well generated 
by the DC magnetic field as illustrated in Fig. 1(b). A 
TT-MW pulse was applied with a randomly chosen phase 
and a normalized amplitude of 7.5 x 10~^ for 500 time 
units to loop 2. This elevates loop 2 to a relatively high 
energy state, while loop 1 is left near the bottom of its 
well. The free evolution of the system is therefore ap- 
proximately given by Eq. (6) with 6d ^ and < B < A 
after the initial tt-MW field on loop 2 is terminated. The 
two coupled oscillators now evolve until a rotation pulse is 
applied to one or both of the oscillators, after which pos- 
sible switching is measured by applying a probe pulse to 
each loop. The rotation pulses were chosen in accordance 
with the experiments with a duration of 200 time units 
and same amplitude as the initial tt-MW pulse. Their 
phases were randomly chosen, but such that X and Y 
pulses were ^ out of phase. The probe pulse was cali- 
brated so that the initial state prior to MW pulses would 
yield vanishing switching probability, the state of loop 2 
immediately following the initial tt-MW pulse would re- 
sult in near certain switching, and the states, where the 
energies of the two loops become equal, would result in 
some intermediate percentage of significance. Thus, we 
chose a half-sine wave with duration 500 time units and 
normalized amplitude 3.6 x 10"'^. The energies of the 
two loops became equal at a time of approximately 900 
units after the termination of the tt-MW pulse, and the 
probe pulse was centered near this value. The {I),X,Y 
rotation pulses were applied immediately prior to the ap- 
plication of the probe pulse. We repeated such simula- 
tions 2000 times with different values of the randomly 
chosen phases in order to obtain a set of switching prob- 
abilities for any given combination of rotation pulses. We 
then generated a density matrix from the simultaneous 
switching probabilities, as outlined above. A typical ex- 
ample of this is presented in Figure 2, where the real 
and imaginary parts of the matrix are shown separately. 
We note the striking resemblance to the density matrix 
that was previously pubhshed (Fig. 3B of Ref. [14]) based 
on the experimental data. The matrix shows relatively 
small real magnitudes outside the diagonal, and relatively 
small imaginary elements, except for p23, which contain 



Re({ft,}) Im({ft,}) 




FIG. 2: State tomography of "entangled" qubits. Real and 
imaginary components, respectively, of the density matrix. 
We deliberately choose the notation rather than to 
distinguish our results as fundamentally classical rather than 
quantum mechanical. This density matrix should be com- 
pared to Fig. 3B of Ref. [14]. 



the signature of " entanglement" . Thus, simulation data 
generated from the classical RCSJ model yields a density 
matrix that matches in its essential features the density 
matrix produced from experimental observations of the 
coupled qubit system. In order to emphasize the similar- 
ity between our result here and that of Ref. [14], we can 

calculate the "fidelity" Foxp = \/ v*"^ PdataV, where pdata 
is the density matrix developed based on the acquired 
classical data and = (Oj ""^j 0)"^ is the vector cor- 
responding to the ideally expected entangled state [5]. 
With the developed density matrix shown in Figure 2 
we get Fexp = 0.77. This value is very close to the one 
obtained in Ref. [14], and we will now discuss this result. 

The probe pulse was calibrated to give approximately 
the same optimal switching probability for the two loops 
when no rotation pulses were applied. The implication of 
this is that the corresponding Bloch vectors would be po- 
sitioned in the xy-plane of the Bloch sphere when these si- 
multaneous measurements are conducted. The probabil- 
ity Pi for switching a loop can therefore be related to the 
z-coordinate of the unit-Bloch-vector in the simple way 
Pi = i(l -|- z). Correspondingly, if an X or a F rotation 
is applied prior to the measurement, one would have the 
probabilities — 5(1 + sin6') and Py^ = ^(1 — cos6'), 
where 6 is the angle (relative to the a;-axis) of the Bloch- 
vector in the xy-plane prior to the application of a ^ 
rotation pulse. Notice that the probability of switch- 
ing in the classical system is given only by the sampling 
of the outcomes of random MW-phases (at the extreme 
low temperature limit), whereas the quantum mechani- 
cal picture sees the switching probability as only due to 
tunneling. In averaging over all random phases 9, one 
can get a good sense of what to expect from the density 
matrix elements. For example the particular element p23 
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corresponding to a random-phase averaged ensemble of 
two Bloch-vcctors, which are T — 6i ^ 62 out of phase 
in the xy-plane, we have P23 — ^exp(ir). This explains 
the information that the density matrix provides. The 
important element P23 is a measure of the angle F be- 
tween the two simultaneously detected Bloch-vectors in 
the xy-plane - or, equivalently, the phase difference be- 
tween the Josephson phase oscillations of the two loops. 
If we produce an clement p23 ~ i f (which is what is ob- 
tained), then this arises from two vectors that are nearly 
orthogonal in the xy-plane - or, equivalently, that the os- 
cillations of the two loops are ^ out of phase. However, 
the formalism does not explain why this is, and the result 
may have nothing to do with quantum mechanics, as we 
have seen here. 

The fact that the classical RCSJ model produces this 
particular result is not a coincidence. We concluded 
above that the tt-MW initiation of the system provides 
for the parameters < B < A and ~ 0. If we ob- 
serve the system at the time tm when the two energies 
coincide {codtm = f)) we can sec from equations (7) and 
(8) that r = 6»i - 6*2 = tan-i ^ - tan"! which 
is r = -f for A = B, and r > -| for A > B > 0. 
Thus, the two coupled oscillators are naturally orthogo- 
nal in phase when their energy contents are similar, and 
the resulting density matrix, obtained from the switching 
measurements after rotations, therefore reflects this fact. 

We have demonstrated that the classical RCSJ model 
of two coupled qubits exhibits the same signatures as ob- 
served experimentally when conducting state tomogra- 
phy and examining the resulting density matrix. These 



signatures have previously been used to argue for an 
unambiguous demonstration of quantum entanglement. 
However, given that our data are definitely classical in 
origin, we submit that such definite conclusions cannot 
be made. Instead, the density matrix, and the element 
P23 in particular, provides information only about the 
phase difference between the Josephson oscillations of the 
two siiperconducting loops without providing the under- 
lying reason for this phase difference. Assuming that 
the system is governed by quantum mechanics this may 
be interpreted as a signature of quantum entanglement. 
However, one does not need to assume this, and our clas- 
sical analysis of the system has revealed that exactly the 
same signature in the density matrix appears as a result 
of the weak coupling. This leads to an ambiguous in- 
terpretation of the observed phenomena. We therefore 
conclude that unambiguous demonstrations of quantum 
behavior in this class of systems must be analyzed not 
only in light of the quantum mechanical model and its 
expectations, but also in light of what the corresponding 
classical model can predict and explain. 



Acknowledgments 

We are grateful for very useful discussions with 
Massimo Bianchi, Maria Gabriella Castellano, Mogens 
R. Samuelsen, and Stuart A. Trugman. One of us 
(J.A.B.) received financial support from the Natural Sci- 
ences and Engineering Research Council of Canada. 



[1] J. Bardeen, Phys. Rev. Lett. 42, 1498 (1979); R. E. 
Thorne, J. R. Tucker, and J. Bardeen, Phys. Rev. Lett. 
58, 828 (1987). 

[2] L. Sneddon, M. C. Cross, and D. S. Fisher, Phys. Rev. 
Lett. 49, 292 (1982); R. A. Klcmm and J. R. Schricffcr, 
Phys. Rev. Lett. 51, 47 (1983). 

[3] S. Pancharatnam, Proceedings of Indian Academic of Sci- 
ence, 44, A, 247 (1956). M. V. Berry, Proceedings of the 
Royal Society of London, A, 392, 45 (1984). 

[4] R. Y. Chiao, Phys. Rev. Lett. 57, 933 (1986); A. Tomita 
and R. Y. Chiao, Phys. Rev. Lett. 57, 937 (1986) 

[5] F. D. Haldane , Phys. Rev. Lett 59, 1788 (1987) 

[6] J. Segort, Phys. Rev. A 36, 10 (1987) 

[7] T. Van Duzer and C.W. Turner, Principles of Supercon- 
ductive Devices and Circuits (Elsevier, New York,1981), 
Chapter 5. 

[8] P. W. Anderson, in Lectures on the Many Body Problem, 
Ed. by E. R. Caianiello (Academic Press NY 1964), Vol. 
2, p. 132. 

[9] A. O. Caldeira and A. J. Lcggctt, Phys. Rev. Lett. 46, 
211 (1981); A. J. Legget and A. Garg, Phys. Rev. Lett. 
54, 857 (1985). 

[10] I. Affleck, Phys. Rev. Lett. 46, 388 (1981). 

[11] J.M. Martinis, M.H. Devoret, and J. Clarke, Phys. Rev. 
Lett. 55, 1543 (1985). 



[12] See, e.g., J.M. Martinis, S. Nam, J. Aumcntado, and 
C. Urbina, Phys. Rev. Lett. 89, 117901 (2002); R.W. 
Simmonds, K.M. Lang, D.A. Hite, S. Nam, D.P. Pappas, 
and J.M. Martinis, Phys. Rev. Lett. 93, 077003 (2004); J. 
Claudon, F. Balcstro, F.W.J. Hckking, and O. Buisson, 
Phys. Rev. Lett. 93 (2004); D. Vion, A. Aassimc, A. 
Cottet, P. Joyez, H. Pothier, C. Urbina, D. Esteve, and 
M.H. Devoret, Forschr. Phys. 51, 462 (2003). 

[13] N. Gr0nbech-Jensen et al, Phys. Rev. Lett. 93, 107002 
(2004); N. Gr0nbech- Jensen and M. Cirillo, Piiys. Rev. 
Lett., 95 067001 (2005); J. E. Marchese, M. Cirillo, and 
N. Gr0nbech- Jensen, Phys. Rev. B 73, 174507 (2006); N. 
Gr0nbech- Jensen et al., "Anomalous thermal escape in 
Josephson systems perturbed by microwaves" , in Quan- 
tum Computing: Solid State Systems eds. B. Ruggiero, 
P. Delsing, C. Granata, Y. Paskin, and P. Silvestrini 
(Springer, New York, 2006), pp. 111-119; J. E. March- 
ese, M. Cirillo, and N. Gr0nbech- Jensen, Eur. Phys. J. 
147, 333-342 (2007). 

[14] M. Steffen, M. Ansmann, R.C. Bialczak, N. Katz, E 
Lucero, R. McDermott, M. Neeley, E.M. Weig, A.N. Cle- 
land, J.M. Martinis, Science 313, 1423 (2006). 

[15] J.A. Blackburn, J.E. Marchese, M. Cirillo, and N. 
Gr0nbech- Jensen, Phys. Rev. B 79, 054516 (2009). 



